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Abstract. We discuss Preheating after an inflationary stage driven by the Standard Model (SM) Higgs field non-minimally 
coupled to gravity. We find that Preheating is driven by a complex process in which perturbative and non-perturbative effects 
occur simultaneously. The Higgs field, initially an oscillating coherent condensate, produces non-perturbatively W and Z 
gauge fields. These decay very rapidly into fermions, thus preventing gauge bosons to accumulate and, consequently, blocking 
the usual parametric resonance. The energy transferred into the fermionic species is, nevertheless, not enough to reheat the 
Universe, and resonant effects are eventually developed. Soon after resonance becomes effective, also backreaction from the 
gauge bosons into the Higgs condensate becomes relevant. We have determined the time evolution of the energy distribution 
among the remnant Higgs condensate and the non-thermal distribution of the SM fermions and gauge fields, until the moment 
in which backreaction becomes important. Beyond backreaction our approximations break down and numerical simulations 
and theoretical considerations beyond this work are required, in order to study the evolution of the system until thermalization. 
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INTRODUCTION 

If the Universe went through an inflationary expansion in an early stage of its evolution, it must also have undergone a 
Reheating period afterwards, during which (almost) all the matter of the Universe was created. A simple consequence 
of Inflation is that the number density of any particle species exponentially dies away with the inflationary expansion. 
Therefore, at the end of inflation, the Universe has no particles at all. Only the homogeneous energy density responsible 
for inflation (usually the potential energy of a scalar field, the inflaton) has survived. So, what? 

The hot Big Bang (hBB) theory describes the early Universe as an expanding space filled up with relativistic species 
in thermal equilibrium. Thus, we are left on one hand, with an empty universe at the end of inflation and, on the other 
hand, with a universe full of relativistic particles according to the hBB. So, how can those two periods be matched? 
Obviously, the energy density responsible for inflation - the potential energy of the inflaton - had to get converted 
somehow into particles. The precise epoch in which the primordial inflationary energy was converted into (almost) all 
the particles of the Universe, is known as Reheating. 

How can particles be created out of the energy of one field (the inflaton)? The idea is that by coupling the inflaton 
to other fields of matter, then the potential energy of the inflaton gets converted into quanta of those fields. A 
specific scenario of Reheating will consist of a specific model of particle physics, specifying a particular form of 
the inflaton couplings to other matter fields. Depending on the model, there will be specific mechanisms of particle 
production taking place. For instance, the main mechanism of production might be perturbative decays OJ], parametric 
resonance 10], tachyonic instabilities JH, etc. Whatever the mechanism of particle production, the created particles 
interact among themselves and eventually reach a thermal equilibrium state with a common temperature. That is 
the Reheating Temperature, 7rh, which represents the energy scale at the end of Reheating and determines the first 
moment in which the description of the evolution of the Universe can finally be put within the hBB theory framework. 

In practice, in order to analyze Reheating, we will study the dynamics of a system described by a potential 

V = Vm (X ) + VintiX , fc, , ...) , ( 1 ) 

where Vj n f is the inflationary potential and Vj nt represents the interaction terms between the inflaton % and the rest 
of matter fields of the model: scalar fields ({0 a }), fermions ({y//,}), vector fields ({A^}), etc. In the last few years, 
simplified models based only on scalar fields or, at most, including (non-gauge) fermionic species, have been the main 
target of study. For example, a simple model in which the inflaton % is coupled to another scalar field (j) and to some 
fermionic field y/, can be described by 

V = V M + V int = V M ( X )+g 2 X 2 <l> 2 + fxW, (2) 



where g 2 and y 2 are dimensionless couplings. Unfortunately, even dealing only with two or three fields, the aspects of 
Reheating are extremely complex 01j Hi- Consequently, most phenomena of particle creation have been discovered 
in simplified scenarios like the one described by (12), where the Standard Model (SM) particles or Dark Matter 
(DM) candidates are simply absent. Models incorporating a gauge principle in some sector of scalar fields (ignoring 
fermions), have also been considered. For instance, Hybrid Inflation models where the inflaton % is a singlet of the 
SM, and the symmetry breaking field <t> coupled to the inflaton, is the SM Higgs doublet. Such models have been 
studied [5] through lattice numerical integration of the classical equations of motion obtained from the Lagrangian 

- j£f = S^* 1 ^* + X - v 2 ) 2 + g 2 * 2 4> + 4> + ^Ti{F^F^ v } , (3) 

with g,X dimensionless couplings, F^ v the non-abelian field strength and 3^ = —ieA^ the gauge derivative, 
coupling <t> with the gauge fields with strength e. 

A realistic scenario of Reheating should account for (almost) all the matter of the Universe. Therefore, it should 
really incorporate a complete gauge theory with all kind of fermions and bosons, since our actual understanding of 
particle physics relies on the Standard Model (SM), based on the SU(3) x SI/ (2) x 1/(1) group involving scalar, spinor 
and gauge fields. Of course, the SM is known to be incomplete and there are several extensions of which one can think 
about. But whatever new ingredients might be added, this does not change the fact that SM particles had to be produced 
during Reheating, since the posterior evolution of the Universe cannot be understood if those particles were not already 
present in an early epoch. Let's then discuss a potential way to realize Reheating within the SM. 



THE STANDARD MODEL HIGGS NON-MINIMALLY COUPLED TO GRAVITY 

Supposse we consider the SM Higgs field non-minimally coupled to gravity, such that the Universe is described by 

S = S SM + S HG = S SM + 1 d 4 x l -M 2 R + % , (4) 

where Ssm is the SM action, M„ w 2.43 • 10 18 GeV is the reduced Planck mass, R is the Ricci scalar, and t, represents 
the strength of the Higgs-Gravity (dimensionless) coupling. In the unitary gauge, the Higgs field can be represented as 
<i> = h/ y/2. The Higgs-gravity sector in the Jordan (J) frame then takes the form 



f(h)R-^ v d^hd v h-U(h) 

r- , Pu2\ /o TTfu\ _ ^ (u2 ..2\ 2 



(5) 



f(h) = (M 2 + ^h 2 )/2, U(h) = -(h 2 -v 2 y , (6) 

where U (h) is the usual Electro Weak Spontaneous Symmetry Breaking (EW SSB) potential of the SM Higgs, with 
vacuum expectation value (vev) v = 246 GeV and self-coupling A ~ ff{Q.\). 

As first discussed in [6], if we require this model to be responsible for inflation, the parameters | and A must be 
related as % ~ 10 5 VA. The reason is that while the amplitude of Cosmic Microwave Background (CMB) anisotropies 
fixes the self-coupling of a quartic potential to be A ~ ^(10~ 13 ) CD, the addition of the non-minimal coupling changes 
that condition to A ~ 10~ 10 <i; 2 [6j]. Since we are thinking of the SM Higgs, for which A ~ &(Q.\), we are then forced 
to require | ~ ^(lO 4 ). For more details see @], or the more recent discussions JH]. For a critical viewpoint see H. 
Here we just want to stress that a large non-minimal coupling is a fundamental ingredient if the SM Higgs is to be 
responsible for inflation. Thus, assuming that correctly describes the early Universe at the required energy scale 
for inflation (see next), we will only focus here on the details of Reheating lllCT . lfllil . just after the inflationary period. 

Analyzing Reheating with a (significant) non-minimal coupling in the Jordan frame can be challenging, so let us first 
perform a conformal transformation to the metric, g^ v — > g^ v = Q. 2 gn V . By imposing the condition f(h)/Q. 2 = Mp/2, 
we can this way obtain the usual Einstein-Hilbert gravitational term. From such a condition one then finds the relation 



The conformal transformation leads to a non-minimal kinetic term for the Higgs field Hill , which nevertheless can be 
reduced to a canonical one by redefining the Higgs as 
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(8) 



Thus, considering the transformations (jT]) and ||S), the total action (ignoring the gauge interactions) in the conformally 
transformed frame - the Einstein frame -, is simply 
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(9) 



(10) 



with U (h) given by ([SJ. To find the explicit form of the potential (TlOb in terms of the new variable r, we must find h 
as a function of %, by integrating Eq. ([§). In our case of a large coupling, | 3> 1, it can be shown 111 111 that 



a z (h)^e aKX , (11) 
with a = y/2/3 and K = 1 /Mp. From ( fTTT i. one finds that the Higgs potential ( fTOb in the Einstein frame looks like 
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where we used the fact that v <C M p . During Inflation % ^> M p / a and therefore the inflationary energy density is 
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GUT: 



(13) 



where Mqut = 10 16 GeV is the GUT scale. Thus, we already know the initial conditions for Reheating in this model. 
The Universe is filled with a homogeneous condensate of the Higgs field, whose energy density is around the GUT 
scale. During Reheating, the Higgs condensate will oscillate around the bottom of its potential. Therefore, in order to 
study the production of particles in the presence of a dynamic Higgs condensate, we first have to derive the couplings 
between the Higgs and the rest of the SM particles, in the Einstein frame. 

In the Standard Model, the masses of the gauge bosons and of the fermions are given by the Higgs mechanism, after 
the Higgs acquires a constant vev. In our case the Higgs field will evolve in time, h = h{%{t)), so the effective masses 
of the fermions and of the gauge bosons will then be changing in time as well, like 



m w = m z cos Gw = ~g2Kx{t)) , m f = -yfh{x{t)) , 



(14) 



where 9w is the Weinberg angle 9w = tan - {gi/gi), an d y/> g\ an d g2 are the Yukawa and the C/(l)y and SU(2)i 
couplings, respectively. The pieces of the SM action of interest for us are: the Spontaneous Symmetry Breaking sector, 
responsible for the W and Z gauge bosons masses, the Charged and Neutral Currents, coupling the SM fermions to 
gauge bosons through the 7r, currents, and the Yukawa sector, coupling the SM fermions with the Higgs, 
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By redefining the fields and masses with a specific conformal weight as 



the three pieces of the SM action <fl~5T > preserve their functional form in the Einstein frame. In other words, the 
Lagrangian describing the gauge interactions between the SM fields, can be written in the Einstein frame identically 
as Lagrangian (15[ in the Jordan frame, as long as the redefinitions dT8T>.(fT9b are considered. A key point here is that, 
not only the form of the interaction terms is known from the gauge principle, but also the strength of the couplings is 
known from the high energy particle physics experiments. Most of the work in the literature on Reheating had only 
been focused on models encoding the different mechanisms that could play a role in the process, with the strength 
of the couplings set essentially by hand. However, in the model under discussion, the relative importance of each 
mechanism of particle creation can be exactly asserted, since for the first time we are studying Reheating in a scenario 
in which the field content, the form of the interactions and the strength of the couplings, are all fixed. From this point 
of view, we can say that Reheating in this Higgs-driven inflation scenario, is a realistic scenario of Reheating. 



REHEATING: TRANSFERRING THE ENERGY TO THE SM PARTICLES 

In order to analyze Reheating, let us first find the dynamics of the Higgs condensate. Soon after inflation ends, the 
effective Higgs potential ([T2l can be approximated by a simple quadratic potential around the minimum, 

V(x) = \m 2 x 2 +AV(x), M=^^~10- 5 M P , (20) 

where AV are some corrections which become negligible during Reheating (see next). Thus, at the end of Inflation, the 
Higgs condensate is formed by Higgs quanta of momentum k = (corresponding to a homogeneous field) and mass 
M. This might seem surprising to the reader familiar with the SM Electro Weak (EW) Spontaneous Symmetry Breaking 
(SSB) process, since there the mass of the Higgs excitations around the true vacuum, is known to be \fXvj V2 <C M. 
However, let us remind the that the EW SSB process, as indicated by its name, takes place at the EW scale £ew ~ 10 2 
GeV, whereas here the energy scale is V^ 4 - E Gm ~ 10 16 GeV. We are considering a very large field approximation of 
the SM Higgs SSB potential, \% | ^> v, as induced through the non-minimal coupling to gravity and the requirement of 
inflation. At the inflationary scale, the system does not see the static Higgs vev v, and only the dynamical amplitude (l22l i 
of X will matter, being % v. 

Considering that the Universe expands with scale factor a t q , then the equation of the homogeneous Higgs is 

X + 3Hx + V'{x)=0 => t 2 x + 3qtx + t 2 M 2 [l + 8M 2 (x)]x = 0, (21) 

with 8M 2 (x) describing the non-linear corrections due to the Higgs' self-interactions. In particular, expanding (fT~2T > 
around the minimum, 8M 2 = —f}\x\ + C^ 2 + ^(Z 3 ) j with P = kMp/V6% 2 and y = 7 A /21E, 2 . Let us then consider 
than those corrections are indeed negligible from the beginning of Reheating, i.e. |5M 2 (j)| <C 1. Then we'll see if 
this can be justified a posteriori. Neglecting such terms, the solution to ( |2TI ) is x(t) x (Mt)~ v J v (Mt), with J v (x) the 
Bessel function of order v = (3q — l)/2. Making use of the large argument expansion (Mt ^> 1) of Bessel functions, 
we then find that the amplitude of the Higgs condensate evolves as 

£(f)«X(f)sin(Mf), X(t)oc(Mt)~% (22) 

Thus, the resulting dynamics of the homogeneous Higgs condensate correspond to an oscillatory field with (angular) 
frequency M and decaying amplitude X(t). The energy and pressure densities associated to the x field, as described 
by ( l22l . can be obtained after averaging over several oscillations, as 

2Px ~ (X 2 +M 2 X 2 ) ~ M 2 X 2 [(cos 2 (Mt)) + (sin 2 (Mt))]=M 2 X 2 , (23) 
2 Px~{x 2 - m2 X 2 ) ~ M 2 X 2 [(cos 2 (Mf))-(sin 2 (Mf))] =0, (24) 

Since the averaged pressure is negligible, then q = d\oga(t )/c/logt is forced to be 2/3, as if a matter-dominated 
background would dictate the expansion rate. Using this, the physical solution can then finally be expressed as 

*( f ) ~ £k sin ( 27rA = ^ sin (^') = X U) sin(TTj) , (25) 

where j = (Mt) /% is the number of times the Higgs crosses around zero and is the number of oscillations = j/2. 
Using now the the (covariant) energy conservation equation, p x = —3(d/a)p x , we can then find that Xe = \/8/3Mp. 



Introducing now d25l l into SM 2 , we find that even just after the first oscillation, |5M 2 | <C 0.122. Thus, from the very 
end of inflation, the Higgs effective potential tends very rapidly to that of a (damped) harmonic oscillator, which 
justifies self-consistently a posteriori the approximation |5M 2 | <C 1 used in the derivation of d25l l. 

Let us now move into the details of the transfer of energy from the dynamical Higgs condensate to the SM particles. 
We will first look at the (perturb ative) decay of the Higgs, and realize that before a single Higgs particle decays, 
more interesting (non-pertubative) phenomena take place. We will then show that, as in the standard case of Reheating 
after chaotic inflation, particles are created when the Higgs condensate oscillates around the minimum of its potential. 
However, the produced particles significantly decay into other species during each oscillation of the Higgs and, as a 
result, a new phenomenological interplay of perturbative and non-perturbative effects has to be taken into account. 



Frustration of the Higgs Perturbative Decay 

A natural Reheating mechanism, would be a perturbative decay process of the Higgs quanta into the SM particles, 
right after Inflation ended. In order to have a perturbative decay of a Higgs particle two conditions must be fulfilled: 

1) The Higgs decay rate T ~ j^M has to be greater than the rate of expansion H 2 = ~ s(^) 2 (^f) ■ 

2) There should be enough phase-space in the final states for the Higgs field to decay, i.e. M > 2my,m^. 

From condition 2), we can deduce ifllll that the Higgs condensate should need to oscillate ~ 10 6 times, before it is 
(statistically) allowed to decay into gauge bosons. And the same applies to top quarks. In the case of the decay into 
bottom and charm quarks, this channel is opened only after a few oscillations of the Higgs, while for the rest of quarks 
and leptons, the decay-channel has sufficient phase space from the very end of inflation. In general, the smaller the 
Yukawa coupling of a given fermion species to the Higgs, the less oscillations the Higgs has to perform before there is 
enough phase-space to decay into such fermion species. However, the smallness of the Yukawa coupling implies also 
a smaller decay rate. Unfortunately, for every species of the SM, when there is phase-space for the Higgs to decay into 
such species, the decay rate does not catch up with the expansion rate. Vice-versa, if the decay rate of a given species 
overtakes the expansion rate, there is no phase-space for the Higgs decay to happen. Therefore, during a large number 
of oscillations, the Higgs is not allowed to decay perturbatively into any of the Standard Model fields. Fortunately, 
before any of such decay channels is opened, many other interesting (non-perturbative effects) start taking place. 



Non-Perturbative Production of Gauge Fields 

Consider for instance the interaction of the Higgs field with the Z gauge bosons (11511811 9b . Around the minimum 
of the Higgs potential, such interaction can be approximated as a tri-linear term \x\Z^.Z^, since the effective mass ( fT~9b 
of the Z gauge boson can be approximated for a small Higgs field amplitude as 

"SIM? 



4<i; cos 2 9w 



For the main part of an oscillation of %, in M and, as a result, the Z gauge boson will oscillate many times during 
each oscillation of the Higgs. Of course, the same applies to the W bosons. Moreover, from ( l26| > and (l22b . we can 
see that during most of the time of Higgs oscillation, the effective masses of the intermediate bosons are changing 
adiabatically, verifying the condition \m\ <C m 2 . However, for values of % close to zero, the adiabatic condition is 
violated, and the passage of % through the minimum of the potential can be interpreted as particle production yy. In 
particular, using the fact that the Higgs velocity around zero can be approximated as w MX(j), the violation of 
adiabaticity can be found in correspondence with a region \%\ <%a 111 111 - with 



V ag 2 M P J \4g 2 % 



where (from now on, unless otherwise stated), g = g2, 82/ cos &w for the W or Z bosons respectively Only outside 
this region, \%\ > % a , the notion of particle makes sense and an adiabatic invariants can be defined. The previous 
regions are indeed very narrow compared to the amplitude of the oscillating Higgs, % a ~ 10~ 2 j x ^X(j). Therefore, the 
particle production taking place in such a narrow field region, happens within a very short interval of time as compared 
to the Higgs period of oscillation T = 2k/M, which can be estimated UU] as At a (j) ~ 2x a /\x\ ~ 10 -2 j^^Af -1 <C T. 
Note that, independently of the species, W or Z bosons, many oscillations (TV ~ 10 3 ) will pass before the fraction of 
time spent in the non-adiabatic zone will increase from a 1% to a 10%, as compared with the period T. 

We will now discuss the non-perturbative creation of particles in the non-adiabatic region. Expanding Eq. d25l ) 
around the j-th zero at time Mti = % j, the evolution equation of the bosonic fluctuations can be approximated (see 111 ill 
for more details) as a Schrodinger-like equation of the type 

- W£ -^-\r\W k = K 2 W k , - Z »-J\x\Z k = K 2 Z k , (28) 

where primes denote derivatives with respect to the rescaled time T = Mt, K = and 

,2 f M pY _8U 



™=™**=-izr\jt) -S;> (29) 

are the usual resonance parameters [2]. Therefore, each time the Higgs crosses zero, we can formally interpret such 
event as analogous to the quantum mechanical scattering of a particle crossing an inverted triangular potential. In 
particular, let T and R = 1 — T be the transmission and reflection probabilities for a single scattering in this triangular 
barrier. The number of particles just after the j-th scattering, n k (j + ), in terms of the previous number of particles 
n k(j ;_ ) just before that scattering, can be written as |01 



n k {j+) = C(xj) + [2C( XJ ) + l]n k (j~) +2cos e jy /c(xj) [C( Xj ) + lWn k {r) [n k (r) + 1] , (30) 

where 9j are some accumulated phases at each scattering (which we will discuss later), and C(xj) is related to the 
transmission probability T k {j) for the j-th scattering as IU2I1 

C(xj) = T-\j) -l = n 2 [Ai (-**) Ai' (-xj) +Bi (-^) Bi' (-xj)] 2 , xj ^ jj^- , (31) 

with fl ; - the scale factor at tj = 7ij/M and Ai(z),Bi(z) the Airy functions. Normalizing the scale factor at the first zero 
crossing as a\ = 1, we can then simply write the evolution of the scale factor as aj = j 2 ^ . Thus, the behavior of xj 
with the number of zero crossings goes as °< j~ 1//3 and from here, the natural range for the momenta of the problem, 
(j), can be found in terms of the q resonant parameters d29| i, as 

Xj = \ => K{j)=q ll3 j lft M. (32) 
Considering the situation in which n k (j~) <C 1, certainly true in the first scatterings, then 

bn k (t)~C{xj) = T^{i)-l, (33) 

where we have retained only the first term of Eq. (l30l . This corresponds to the spontaneous particle creation of W 
and Z bosons each time the Higgs crosses around zero. In particular, the total number of produced particles of a given 
species (and polarization), just after exiting the non-adiabatic region around the j-th zero-crossing, can be obtained as 

f juilr^t., >i _ 1 a nth 



A»(; + ) = —=-3 / dkk 1 C( Xj ) = ^ J>M , (34) 
zn^ajJo 2.J 

with .y = n~ l [q C(x)x 2 dx » 0.0046 and q the resonant parameters given by Eq. ( |29] >. Thus, the only difference 
between the number of W and Z bosons produced is simply encoded in the different resonance parameter, qw « g\ and 
qz°< g 2 / cos 2 %. 



1 All the couplings are renormalized at the inflationary energy scale, ~ Mqut, so here g\~g\~ 0.30 and sin 2 9y/ = cos 2 9w 0.7. 



As pointed out first in [2], one would expect that during the first oscillations of the inflaton (the Higgs in our case), 
the first term of (f30b would dominate over the others, and the particle production should be driven by spontaneous 
creation at the bottom of the potential, as described by d34l >. Eventually, after some oscillations, a significant number 
of particles would have been created and, consequently, the terms proportional to in ( f3Qb should dominate. Those 
terms, as opposed to d33l l, induce a stimulated growth of particles in an explosive manner, known as parametric 
resonance. The accumulation of particles is thus expected to enhance the rate of production, entering into a regime 
of exponential growth. At least, that is the usual picture expected for Reheating yfl in a scenario in which the inflaton 
oscillates around its potential. However, in our scenario under discussion, a new phenomenon is going to occur. 



Perturbative Decays of Gauge Bosons 

In the Reheating scenario under analysis we have access to the form and strength of the couplings between all 
particles. Therefore, we can compute the total decay widths of the W and Z bosons into any pair of fermions. As 
opposed to the conventional result found in particle physics books, now the gauge bosons' masses are time-dependent 
functions of the dynamical vev of the Higgs and, as a consequence, the decay widths will vary in time during Reheating. 
The functional form of the decay rate in the Einstein frame will preserve the same form as in the Jordan frame, being 
only changed through the conformally transformed masses, see Eq. ( fT9] >. After some calculus, one finds 

£ _ 3g 2 2 m w _ 3g 3 2 M p ( | |\i/2 rE _ 2Lips rE 

Fw± ~~l6r~32^P-V' e J ' rz_ 3co S 3^ rw± (35) 

where Lips = | — -U- sin 2 $w + t sin 4 @w denotes the Lorentz invariant phase-space factors. Note that after the passage 
of x through the bottom of the potential, the gauge bosons' masses grow as the Higgs climbs up the potential and, 
as a result, the decay width also increases significantly. The W and Z bosons then tend to decay into fermions in a 
time At ~ 1 /(r^ z ) ; -, where (•)_,■ represents a time average between the j- and the (J + l)-th Higgs zero crossings. 

The typical time of decay turns out to be At ~ e j l / 2 M~ 1 , with e = 0.64 or 1.55 for Z and W bosons, respectively. 
Therefore during the first oscillations of the Higgs, the non-perturbatively produced gauge bosons significantly decay 
within a semiperiod T/2= nM . As the amplitude of the Higgs field decreases with time, the decay width of the 
gauge bosons (|35T > also decreases and becomes less and less significant. However, as we will explain next, simply the 
fact that bosons decay during the first oscillations of the Higgs, will have important consequences. 



Combined Preheating: Mixing Perturbative and Non-Perturbative effects 

The total number density of gauge bosons n(j + ) present just after the j-th crossing will decay exponentially fast 
until the next crossing, due to the perturbative decay into fermions. Therefore the total number density just previous 
to the (jf + l)-th zero crossing, n((j + 1)~), is given by 

«(C/ + ID = n(j + y J ' J J +lm = n(j + )e-Wjl . (36) 



On the other hand, in the large occupation limit «a ^> 1, we can neglect the first term in Eq. ( l3Qb - We then find that 
the spectral number densities of the gauge bosons just after and previous to the j-th scattering, and nk(j~), 

respectively, become proportional to each other as 

n k (j+) a ((2C(Xj) + 1) - 2cosdjyJc(xj)(C(xj) + 1)) n k (D = n k (j-)e 2 *^ , (37) 

where C(xj) was defined in OTb and /^(j) is the Floquet index. The {9j} are some accumulated phases at the j-th 
scattering, which can indeed play a very important role, since they can enhance (cos 6j < 0) or decrease (cos 6j > 0) 
the effect of particle production at each crossing. In particular, we can estimate [11] them as 

A6j = dfjk- + m-a *jj&-F(j) ~ <?(10 2 )r "\ (38) 

f «4f"" i '( , -'"" l " wl )^ a4 ''"" 1 ' <39 » 



When A9j = — 9j ^> it, the effect of resonance will be chaotic, such that then the phases essentially are random 
at each scattering. From fl38l l, we conclude that such stochastic behavior occurs for the first <^(10 3 ) oscillations. Thus, 
for the first thousand oscillations the successive scatterings are incoherent, what allow us to define for each scattering, 
an average Floquet index obtained as 2%^ = Jq Z /J./ C (9)d9, see lUlll for more details. 

On the other hand, the decay widths (l35l l of the gauge bosons, averaged over the Higgs oscillations, are given by 

/r \ ( 82 V ^fLips j,, .n _ 2 7z F / .s /r \ 3cos 3 tv _2jw 

{rz - a,l) i = ) J^/I F{J) = — F{]) ' {TW - U) i = ~2L^- {Tz -» u) J = — F{]) ' (40) 
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where F(j) is defined in ( f39l > and the constants Jz,Yw are just numerical factors depending of the parameters of the 
model and the decaying species. The decay of the vector bosons occurs precisely between two successive Higgs zero- 
crossings. Thus, taking into account Eqs. d37l ) and d36l l, the number of gauge bosons just after the (j + 1 )-th scattering, 
in terms of the number just after the previous one, can then be expressed as 

«*((./' + 1) + ) =B Jt (C/'+ir)« 2,t * C/+1) = n t 0" + )e~ rF( V*' J * (; ' +1) , (42) 

This formula ( l42l clearly shows how the two effects, the non-perturbative resonance production (via e 2z ^) and the 
perturbative decays (through e~^ F ^), are combined together. This combination indeed generates new phenomenology 
and to emphasize the difference from the usual parametric resonance, we will call it Combined Preheating. Apply- 
ing (l42l recursively, we could estimate the occupation number (for each species and polarization) just after the y'-th 
scattering, n/ ( (j + ). That was indeed the procedure we adopted in 111 ill . However here, to be more accurate, we will 
rather use the recursive iteration of (f43t . The reason is that we obtained (l42l by neglecting the first term of d30i l. which 
is a reasonable approximation when parametric resonance takes place very soon after the end of inflation (as implicitly 
assumed in chaotic models). However, in the model we are discussing, the gauge bosons decay significantly. Thus, 
they are intially not sufficiently accumulated and, consequently, resonance is delayed. In particular, considering (l30b 
and d36i > all together, one obtains the phase averaged relation 

Q +»*((; + 1) + )) =(l+2C(xj))(±+n k (j + )e-^j)\ , (43 ) 

where we simply dropped out the term involving cosQy, after averaging Eq. (l30b over the random phases {Oj}. We 
clearly see that the perturbative decays e~^ F ^' tend to decrease the rate of production of bosons, while the factors 
(1 + 2C(xj)) > 1 [or equivalently the factors e 2rt ^ k in (l42l l tend to (resonantly) enhance it, due to the accumulation 
of previously produced bosons. Initially, the perturbative decays prevent the resonance to be effective. However, after 
a certain number of oscillations, the resonant effect may will finally be developed, since the perturbative decays 
become less and less important as times goes by. That is precisely the novelty in Combined Preheating, a simultaneous 
competition between perturbative decays and non-perturbative production of particles, affecting each other recursively. 
Depending on the couplings, the energy density transferred into the fermionic species might dominate over the Higgs 
condensate before parametric resonance becomes relevant or, perhaps, the energy transferred into the fermions is 
negligible and the perturbative decays simply delay the development of parametric resonance. Since in this model we 
know the strength of the couplings, we will univocally determine which of the two previous options is the correct one. 
However, note that Combined Preheating might be a common ingredient (ignored so far) to many Reheating scenarios. 
As soon as the inflaton is coupled to some primary fields and the latter to some secondary fields, a competition between 
the non-perturbative production of the primary particles and a pertubative production of the secondary (from the decay 
of the primary ones) will take place. The rate at which energy is exchanged between the inflaton and the primary and 
secondary fields, could only be determined in each specific scenario if the involved couplings were known. However, 
the couplings of the inflaton to other matter fields are essentially free parameters in most models. Therefore, in the 
light of all this, we plan to explore in the future the space of parameters of other models different than the one under 
discussion, in order to determine the role that Combined Preheating could play. 

Coming back to our model, it can be shown ll ill that both W and Z bosons are non-relativistic, while their fermionic 
decay products are all ultra relativistic. Thus, the mean energy of the fermions (F) produced from the decay of the 
gauge bosons (B) between tt and tj+\, can be estimated as 

E m (j) ee ( (4+m!) 1 / 2 ); w (k F )j w 0.5(m B )j a 0.25(g/^ 2 )F(j)M p , (44) 




FIGURE 1. Evolution of the energy density transferred into the gauge bosons and into the fermions as a function of j, for X = 0.2 
and | = 44700\/X The decay of the homogeneous energy density of the Higgs is also shown. All densities are in units of M 4 . 



Then, the number of fermions produced between the j-th and the ( /' + 1 )-th scatterings and their corresponding energy 
density, are simply given by 

An F (j) = 2x3 x (»z(; + )(l-^ F(/) )+2« w (j+)(l- e -^)) , (45) 



Ap F (J) =2x3 [(1 -e-^U) )nzU+)EF(z]U) +2(1 -e-^j) )nw{j + )EFmU) 



(46) 



where the factor 2x3 takes into account that each gauge boson can have one out of three polarizations and decay into 
two fermions, while the extra factor 2 in front rtw, accounts both for the W + and W~ decays. After j crossings, the 
total energy density transferred into the relativistic fermions and into the non-relativistic bosons, can be expressed as 

PF(j) = i(i/j) &/3 Ap F (i), (47) 

Pb(J) = 3 {n Z (j + ){mz)j+2n w (j + )(m w )j) , (48) 

By iterating recursively (l43l we can find «z(j + ) and nw(j + ) at eacn crossing, and plug them into d46b.(l47]i and (08]). 
This way, we can follow the time evolution of the energy density of the Bosons and of the Fermions. The results 
are summarized in Fig. (Q]). We find that the number of zero-crossings jg for which the perturbative decays stop 
blocking resonance, is jx ~ 70 for the W bosons and jn ~ 300 for the Z bosons. Parametric resonance thus becomes 
important much earlier for W's than for Z's since their decay rate ( 1351 ) differ in a factor Jz/jw ~ 2.4, such that many 
more W bosons survive per half period than Z bosons. Therefore, the Combined Preheating of the W bosons is much 
faster driven into the parametric-resonant like behavior, while the evolution of the Z bosons is much more affected 
by the perturbative decays, completely preventing the development of parametric resonance. After several dozens of 
oscillations, the transfer of energy from the Higgs to the gauge bosons is completely dominated by the channel into the 
W bosons, since they become fully resonant while the Z bosons are still severely affected by their perturbative decay. 

We can estimate the time in which finally the energy of the inflaton would be transferred efficiently to the fermions 
or the bosons. Defining that moment, respectively, like e F (j e s) = p F /p x = l ar >d £ b C/eff) = Pb I P% = 1 > one obtains the 
numbers in Table Q] Note also that the number of oscillations j e ff required for an efficient transfer of energy, depends 
on the parameter A, although the overall order of magnitude does not change appreciably. 

Unfortunately, before reaching the stage in which e F b ~ 1, the backreaction of the produced gauge fields into 
the homogeneous Higgs condensate becomes significant. We don't have space here to give a detail account of 
this, so we will simply summarize our findings in ifllll . There we found that the effective oscillatory frequency of 
the Higgs, once the Higgs - Gauge fields interactions are restored in (f2TT >. is given by an expression of the type 
ft'osc = M 2 [l + 0(gj 3 / 2 /i; i / 2 )(nw(j)/M 3 )]. Thus, since resonance is initially blocked, the second term is very small 
and ©osc w M, as we have implicitly assumed all the time. However, when parametric resonance becomes efficient, % 
grows exponentially fast within few Higgs oscillations, and the term in co 2 sc proportional to nw eventually dominates. 
We can determine the number of Higgs crossings, jbr, such that for j > j^, backreaction of the bosonic fields cannot be 
ignored anymore. The results, as a function of A, are summarized in Table Q] We clearly see that backreaction seems 



TABLE 1. Left: Number of Higgs crossings (depending on X) for efficient Reheating. Right: 
Number of Higgs crossings (depending on X) for efficient Reheating. 
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to become important at a time slightly earlier than that at which we were expecting the Higgs to have transferred 
efficiently its energy to the bosons and fermions. This means that our analytical estimates of these transfers were 
biased, and a careful numerical study of the process is required beyond backreaction. 



FINAL REMARKS 

We have studied Preheating after an inflationary stage driven by the Standard Model (SM) Higgs field non-minimally 
coupled to gravity. We have found that Preheating is driven by a complex process, which we called Combined 
Preheating, in which perturbative and non-perturbative effects occur simultaneously. 

The Higgs field, initially an oscillating coherent condensate, produces non-perturbatively W and Z gauge fields. 
These decay very rapidly into fermions, what prevents gauge bosons to accumulate and, consequently, blocks the 
usual parametric resonance. The energy transferred into the fermionic species is, nevertheless, not enough to reheat 
the Universe, and resonant effects are eventually developed. Soon after resonance becomes effective, also backreaction 
from the gauge bosons into the Higgs condensate becomes relevant. We have determined the time evolution of the 
energy distribution among the remnant Higgs condensate and the non-thermal distribution of the SM fermions and 
gauge fields, until the moment in which backreaction becomes important. Beyond backreaction our approximations 
break down and numerical simulations and theoretical considerations beyond this work are required, in order to study 
the evolution of the system until thermalization. 
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